Phase diagram and quench dynamics of the Cluster-ZF spin chain 
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In recent years, the understanding of quantum phases and 
quantum phases transitions has gained an incredible momen- 
tum. This is due to both strong theoretical advance, stimulated 
by the richness of quantum phases, and experimental advance. 
On the theoretical side, there has been great progress in under- 
standing how quantum critical points affect the finite temper- 
ature regime Q and quantum entanglement| 2 1 . Moreover, we 
are realizing that quantum phases host rich novel quantum or- 
ders and phases of matter |3|. From the experimental side, 
ultra cold atom gases have proven to be the ideal arena to see 
coherent quantum evolution for many body systems (UQ. 

The study of closed quantum systems out of equilibrium is 
important for manifold reasons. To start with, one is inter- 
ested in applications to quantum information, in which deco- 
herence and entanglement dynamics play a fundamental role. 
An important theoretical perspective regards the understand- 
ing of the notion of universality for a system away from equi- 
librium, where the traditional concepts of phase, renormal- 
ization group and fixed point fail. Recently, the study of the 
equilibration of quantum many body systems has given new 
insights in the foundations of statistical mechanics ITtUToII. 

Driving a system out of equilibrium can be accomplished 
in many ways. Most of the efforts have focused on quan- 
tum quenches ifTOl , namely sudden global or local changes 
of the external parameters of the Hamiltonian governing the 
unitary evolution of the closed system. One of the ways of 
understanding the dynamics of a system after a quench is the 
Loschmidt echo, a measure of the partial recurrences with the 
original state as a function of time |[IIl[l2 |. Recently, the time 
behaviour of the Loschmidt echo has been investigated in var- 
ious models, in particular the XY spin chain ri3llT6l. 

In this paper we will consider a one-dimensional model 
that extends the XY spin chain with a three-body cluster term. 
The exact solution becomes available using well-known tech- 
niques and this allows us to study the complete phase diagram. 
We find that a particular critical region has a behavior quite 
different to the one found in the XY model. We then study 



the behavior of the Loschmidt echo after critical quenches for 
two different critical points. We find qualitative differences 
derived from the non trivial nature of the phase space. 

Cluster-XY spin chain. — Cluster states have emerged re- 
cently as a physical system for implementing one-way quan- 
tum computation. In particular, it has been shown that two- 
dimensional cluster states serve as fiducial states for universal 
measurement-based quantum computation I TtHITII . 

The cluster state can be defined using a so-called sta- 
bilizer Hamiltonian. Consider a finite-dimensional lattice 
L composed of N vertices, each vertex containing a two- 
dimensional quantum system (qubit). The cluster state for this 
system can be defined as the unique +1 eigenstate of the sta- 
bilizer operators 



n 



(1) 



where v ^ ji denotes that v is connected to ji and a" are the 
Pauli matrices [22 1 . The stabilizer Hamiltonian is simply 



c 



and the cluster state is defined as its ground state. This 
preparation can also be achieved by preparing all the qubits 
in the |0) state (a^ |0) = |0)) and then performing a con- 
trolled sign operator U = exp(i7r |+) (+| |+) (+|) (where 
cr^ 1+) = 1+)) on every pair of connected vertices ll23ll . The 
model we study incorporates the one-dimensional version of 
the cluster phase competing with the XY model in a transverse 
field. The Hamiltonian is 
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where cr^{a = x^y^z) are the Pauli matrices acting on the site 
n of the lattice and we impose periodic boundary conditions 
(^A^+i — ^?)- Similar models were considered in 



2 



Defining local raising and lowering operators cr^ = ^ (cr^ ± 
io-y), we obtain global canonical anticommutation relations 
by using a Jordan- Wigner transformation 1 1 1 



4 = {U<] 



(3) 



SO that the model is mapped to a quadratic Hamiltonian of 
spinless fermions {cn,c^} = 0, {cn,cj^} = Snm- Note that 
the parity operator Q := Yl^ commutes with the Hamil- 
tonian and can be diagonalized simultaneously with it. Using 
the fact that the system has translational invariance, we may 
perform a Fourier transform 



Ck 



1 ^ 

^ n=l 



N 



(2m + l-g), 



where we decompose the Hilbert space so that Q = (—1)^ 
and m = 0,...,A/' — 1. We can then rewrite the Hamiltonian 
as 

0</c<7r 

up to a constant, where 

Ck = cos(2/c) — {Xx + Xy) cos{k) — /i, (4a) 
5k = sin(2/c) - {Xx - Xy) sin(/c). (4b) 

We diagonalize the Hamiltonian by means of a Bogoliubov 
transformation 



7/e = COs((9/e/2)c/e - Z sin((9/e/2)cL;,, 



(5) 



imposing 6_k = -Ok so that {7fe, Tfc/} = The Hamilto- 
nian becomes 

H = 2 ^ Ak (tIt/c + tI/cT-Zc - l) + const, (6) 

0<fc<7r 

where we defined the en ergy for the so called Bogoliubov 
quasip articles = + and we canceled the un- 

wanted 77 terms by choosing Ck sin 6>/e + 5k cos 6>/e = 0, or 
equivalently 



Ok 



arctan 



(7) 



The ground state has the form of a BCS state in terms of the 
original operators 

l^)= n (cos(^fc/2)+zsin(^fc/2)4cLfc) |0)^, (8) 

0<fc<7r 

where Ck |0)^ = 0, V/c. 

Phase diagram. — At this point, we can draw the phase di- 
agram by finding the regions of quantum criticality where the 
system becomes gapless in the thermodynamics limit, that is. 
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Figure 1. Reduced phase diagrams for (a) A^c = and (b) /i = 0. 
We use the conventions: P (polarized, given by cr^), C (cluster), 
AFM (antiferromagnetic), FM (ferromagnetic). 



A/e = for some k G [— tt, tt) when N ^ oo. First, note that 
trivially ^/^ = for /c = 0, ±7r. In that case, ek vanishes for 



h = ±{Xx^Xy)^l. 

Now, if A: ^ 0, TT, J/e vanishes if cos(A:) = 
relation for e/c, we get the critical region 



(9) 

. Using this 



h — Xy — XxXy 



1, 



-2 < A. 



A^ < 2. 



(10) 



We see an asymmetry between A^^ and Xy . This should be no 
surprise as this is already evident in the cluster term of the 
Hamiltonian. Note also that the union of the critical regions is 
invariant under Xy ^ —Xy,Xx ^ —Xx. We can understand 
this by noting that both local unitary transformations 



'2n 



acting only on even sites and 



^2 • 0-2n-^i ^ <^2n+l5 ^2n+l ^ <^2n+l' ^2n+l ^ ^2n+l 

acting only on odd sites map H{Xx^Xy,h) to 
H{—Xx^—Xy^h). This is a consequence of the Z2 x Z2 
symmetry of the cluster state implemented precisely by Ui 
andf/2 1231. 

One of the interesting features of this model is the existence 
of phases that appear because of the competition between the 
XY and cluster terms. Consider first A^^ = (Fig. 1(a) ). The 
Hamiltonian in this case does not have Ising interactions of 
the type cr^cr^+i- However, two of the regions next to the 
cluster phase can be connected adiabatically to ferromagnetic 
and antiferromagnetic states in the x direction, respectively. 
Something similar happens fork = 0. The Hamiltonian does 
not have a transverse term that tries to polarize all the spins in 
the same z direction, nevertheless this phase is present in the 



reduced phase diagram (Fig. 1 1(b) ). 

Fidelity susceptibility. — The phase diagram can be stud- 
ied [261 by considering the fidelity susceptibility introduced in 
1271 . namely the response of the ground state to small changes 
of the external parameters. Consider a many-body system de- 
scribed by a Hamiltonian 



H{X) = Ho^XHi, 



(11) 
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where A is an external parameter used to control the system. 
There is no loss of generality if we write the Hamiltonian this 
way, especially if the system is large enough so that the critical 
point is well localized. Hi is called the driving Hamiltonian. 
We may now diagonalize the system and obtain both the en- 
ergy spectrum and the eigenstates 



H{\)\n{\))=En{\)H\)). 



(12) 



If we change A to A + (5A and we are away from possible crit- 
ical points, the physics described by the neighboring ground 
states will be similar. However, note that in the thermody- 
namic limit different ground states will become orthogonal, 
as was realized by Anderson in the so-called "orthogonality 
catastrophe" |28|. For finite size systems, we expect that the 
new ground state will remain "close" to the original ground 
state and we may study how fast does the overlap goes to zero. 
In order to quantify this notion we use the fidelity EZl ESI 



^(A,A') :=|(1^(A)|1^(A'))|, 



(13) 



where |f2(A)) represents the ground state of H{\). The re- 
sponse of the fidelity after an infinitesimal change of the ex- 
ternal parameter up to second order reads 



a2 



A + <5A) = 1 



where the fidelity susceptibility I29H3TI is defined by 



(14) 



- {dxn{\)\n{x)) dx^ix)) . (i5) 

If we have more than one external parameter, we may gener- 
alize this result and obtain the so called quantum geometric 
tensor ISTIl 



Tab ■.= {dxM>^)\d>,,^{X)) 

-{d>,^^l{X)\^l{X)){n{X)\dx,^{X)). 



(16) 



In general, this tensor will be complex. Both the real and 
imaginary parts have nice physical interpretations 126]. The 
real part will be an induced Riemannian metric on the mani- 
fold of parameters 



9ah (Tab). 



(17) 



The geodesies with respect to this metric give information 
about the optimal adiabatic path that connects to points inside 
the same quantum phase. Also, the scalar curvature may be 
used to distinguish different phases and characterize the be- 
havior of critical regions |26|. The imaginary part is related to 
the appearance of geometrical phases through the Berry cur- 
vature 

Fab :=9(Ta6) 

= {d,M><)\dxM><)) - {dxM^)\dxM><)) 

= daAb - dbAa, (18) 

where A := d\^D,) is the adiabatic Berry connection lISTl 
|32J. 
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Figure 2. (Color online) Contour plot of the fidelity T{\y , \y+5\y) 
SXy = 0.05, N = 500, for constant h: (a) /i = -1.5, (b) /i - -1 
(c) /i = 0, (d) /i = 1 and (e) /i = 2. (f) Constant A^^ = 



Coming back to the cluster-Ising model, we may compute 
the fidelity between neighboring ground states |^^(A^)) ([5]) for 
different values of the {A^}. If we change the external pa- 
rameters A-^^ x[^\ we can express the "old" ground state 
1](A -^^ ) ^ in terms of the operators that diagonalize the "new" 

(2) 

Hamiltonian H{Xl The form of the wave function remains 
the same 



k>0 

(19) 

where 7^ , 7^ are the fermionic operators that diagonalize the 

(2) 

Hamiltonian -ff(A- ') and we defined 



Xk = OkiX\ 
After a straightforward calculation, we obtain 

Ok - 



(20) 
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Figure 3. (Color online) (a) Overlap between |Q(P1)) and the neigh- 
boring ground states (N = 500, ^ = 1 and /i = 0). (b) Overlap 
Fi(A) between the two particle states of PI and the neighboring 
ground states (N = 500, q = 1 and h = 0). 
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Figure 4. (Color online) (a) Overlap between \ Q{P2)) and the neigh- 
boring ground states (N = 500, q = I and /i = —1). (b) Overlap 
Pi (A) between the two particle states of P2 and the neighboring 
ground states (N — 500, q = 1 and h = —1). 



From this expression, we compute the quantum geometric ten- 
sor EB 

0<k<7T ^ ^ 

using the convention Ai = A^^, A2 = and A3 = h, where 



dOk cos{k)Sk — sm{k)ek 

dOk cos{k)Sk + sin(/c)e/e dOk 



A? 



dh 



Note that Tab may not be analytic when the system becomes 
gapless, i.e., when ^ for some k. Since Tab is a real 
tensor, this system will have a trivial Berry curvature. The 
above expressions for the quantum geometric tensor are non 
trivial. We expect to find a richness of features in their scaling 
behavior [ 31 1, which can be potentially of use in the optimiza- 
tion of quantum adiabatic algorithms |33|. A thorough study 
of the scaling of the geometric tensor in the Cluster-XF model 
is to be found in |34|. 

By computing the fidelity in the Cluster-XF model, we find 
the expected critical lines. We illustrate this in figure [2] plot- 
ting the phase diagram region that we already discussed in the 
context of the exact solution. 

We notice that there is a set of multi-critical points that 
present anomalous behavior. It is known that some multi- 



critical points may behave differently, giving rise to anoma- 
lous dynamical scaling properties and, as a consequence, dif- 
ferent universality classes |[35l|36l. The multi-critical points 
A*^^^ , given by 



h-3 h- 



\/h e 



(21) 



have properties that are not present in other multi-critical 
points, like the point (A^c, Xy^h) = (0, 1, 0) that was studied 
extensively in ll25ll . In fact, the overlap between the ground 

state state corresponding to the critical points 1](A-^^)^ 

and the neighboring ground states, is very large. More- 
over, also the overlap i^i(A-) between such multi-critical 
ground states and the subspace generated by excited pairs 



is quite large: 



0<fe<^ 



(22) 



This overlap region will roughly follow the truncated surface 
given by Eq.([TO). 

We can gain some physical intuition about these multi- 
critical if we rewrite the Hamiltonian as 

+ (h+s)(±^- E(«+i + ^r^r+i) - E • 

This corresponds to the sum of three critical Hamiltonians, 
namely a critical Ising model, a critical cluster-Ising model 
and a critical XX model with a transverse field. Note that for 
\h\ 1, the XX model term dominates. The ground state of 
this Hamiltonian corresponds, after a Jordan- Wigner transfor- 
mation, to a completely empty (full) Fermi sea QJ • This state 
is trivial from the entanglement point of view since it corre- 
sponds to a product state 1 37 1, but it affects the behavior of the 
fidelity susceptibility in its vicinity 1381 [39l . 

For the sake of concreteness we will concentrate on two of 
these points 



PI : 



P2 : 



(Ax, Xy^h) — 



3 1 

2' 2 



,0 



(A„A^,/i) =(-2,0,-1). 



(23) 
(24) 



The overlap with the neighboring ground states is illustrated 
in figure |3] and |4] We see a significant overlap between the 
state at the quantum critical point and the neighboring ground 
states of the antiferromagnetic phase. This phenomenon is due 
to the fact that the perturbation corresponding to the parameter 
A^ is not sufficiently relevant |31 1. The overlap i^i(A-) ([22]) is 
also considered, showing that the most significant overlap is 
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Figure 5. LE starting from the cluster phase to the critical point Xx = 
- 1, Ay = I and /i = with N = 400 starting from (A^ and h kept 
fixed) (a) Xx -1.3, ^ = 1, (b) A^ -1.3, ^ = 0, (c) A^^ -1, 
^= l,(d) A,. = -1,^ = 0. 



either with the ground state or just a few pairs of excitations. 
Note that Fi is symmetric, so we also get the overlap between 
the subspace of a pair of excitations of the neighboring states 
and the critical point ground state. 

The other critical points that do not belong to these criti- 
cal lines present a behavior similar to the one obtained for the 
XY model in previous studies ll27l . In those cases, the over- 
lap with all the neighboring ground states decays very fast 
even for finite systems. In the present model, we expect this 
behavior for asymptotically large values of all the couplings 
Xx, Xy and h since the cluster interaction in the Hamiltonian 
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Figure 6. LE for the quenched cluster-Ising model to the critical 
point Aa, = - 1, Ay = I and /i = with N = 400 starting from (a) 



Ay = 0.7, ^ = 1 (Xx and h kept fixed, polarized), (b) Aa; = 
q = 1 (Xy and h kept fixed, ferromagnetic). 



-1.7, 



^ becomes negligible in comparison and we obtain the usual 
XY model. This fact implies that the two planes given by ^ 
correspond asymptotically to the Ising critical lines and the 
hyperbolic surface ([TO]) correspond asymptotically to the XX 
critical line. Note, however, that the universality class may 
be affected as we get close to the multicritical region we have 
been considering ||36|. 

Quantum quenches and Loschmidt echo. — At this point, we 
are ready to study the dynamics of the system after a quantum 
quench. In order to quantify this we use the Loschmidt echo 
(LE). This quantity is widely used in many-body physics, in 
particular in the field of quantum chaos Iillttl4il . Suppose we 
want to compare the dynamics under the Hamiltonians Hi and 
H2 (possibly time dependent) imposing the same initial con- 
ditions |?/^(t = 0)) = IV^o)- In that case, we define the LE as 



(25) 



where Ua(t) = Texp(— z Ha(t')dt') and f denotes time 
ordering. Note that C{t = 0) = 1. In this paper, we will 
limit ourselves to ground states of one of the Hamiltonians, 
so that one of the unitaries in ( [25] ) acts trivially. This can 
be interpreted operationally as preparing the system in the 
ground state of the Hamiltonian with parameters x[^\ sud- 
denly switching A-^^ A-^\ and letting the system evolve 
with the new Hamiltonian. The LE reads 



(26) 



In this sense, the LE is a dynamical version of the ground 
state fidelity. High values of the LE mean that the system 
is approaching the initial state. Typically, the LE will decay 



6 



exponentially at first and then start oscillating around a well- 
defined value |[l3l[T4l. If the systems is finite, we expect the 
time evolution to be quasiperiodic, driving the system arbitrar- 
ily close to the initial state for long enough times. The system 
will experience revivals, i.e., times where the value of the LE 
is greater than two standard deviations from the average. The 
structure of these revivals may be greatly affected by critical- 
ity |15|. The time evolution in the Cluster-XF model after a 
quantum quench is given by 



m)) 



n (cos(f ) + ze--^^ sin(f )7l7l.) \^{>^P)) ■ 



0</c<7r 



and the LE is thus 



(27) 



c{t) = \m)\m)? 

= Yl (l-sin2(xfe)sin2(2tAfc)). (28) 



0<k<7r 



In the following, we show the detailed analysis of the LE for 
different types of quenches. 



Consider the critical point h = 0, 
critical Hamiltonian is 
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2 = 1 
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As we see in figure 1(b) this point lies between three different 
phases. If we increase A^; (A^) we will be in the cluster (po- 
larized) phase. Decreasing \y or A^^ results in a ferromagnetic 
state in the x direction. 

For this critical point the behavior of the LE depends on 
which phase the system is prepared in. For the cluster state 
(Fig. [5| the numerical simulations show that it will oscillate 
strongly away from the mean value. We also get a noticeable 
spreading of the revivals. If we start further from the critical 
point (i.e. a stronger quench), the lines get closer, but we keep 
the insensitivity to the parity of the ground state (Q = (—1)^). 
Note that the stronger the quench, the sooner the revivals will 



happen. If we start from the polarized phase (Fig. 6(a)), the 
behavior of the LE is somewhat different. Once again, the 
oscillations make the revival structure insensitive to the parity 
sector. 

On the other hand, starting from the ferromagnetic phase 
changes the LE completely (Fig. |6(b)|). The numerical simu- 



lations show that it will oscillate randomly around a relatively 
high mean value. There is no outstanding structure for the 
revivals as we had in the previous quenches. Increasing the 
size of the quench will give us basically the same result, only 
decreasing the mean value. This is consistent with the signif- 
icant overlap of the ground state of this critical point with the 
neighboring ferromagnetic ground states (Fig. |3(a)| ). 

Consider now /i = 0, Xx = 0, Xy = 1. For these values, 
the critical Hamiltonian is 
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Figure 7. LE for the quenched cluster-Ising model to the critical 
point Xx = 0, Ay = 1 and /i = with N = 400 starting from (A^ 
and h kept fixed) (a) A^ = 0.8, q^l (cluster), (b) A^ = 0.8, ^ = 
(cluster), (c) Xy = 1.2, ^ = 1 (antiferromagnetic), (d) A^ = 1.2, 
^ = (antiferromagnetic). 



This critical point lies on the interface of four different 



phases (Fig. 1(b)). If A^ is increased, we will be in a region 
that can be connected to an antiferromagnet in the y direction. 
If A^ is decreased, we will be in the cluster region. A positive 
Xx will turn the system into an antiferromagnet in the x direc- 
tion, while a small negative A^^ will put it in a region that can 
be connected to a separable state polarized in the z direction. 

The behavior of the LE after a quench of Xy (Fig. [7]) is sim- 
ilar to the one obtained for critical quenches in the XY model 
|[T6l . It features the same sensitivity to the parity that cancels 



7 




100 200 300 400 500 
t 

(a) 




100 200 300 400 500 
t 

(b) 



Figure 8. LE for the quenched cluster-Ising model to the critical point 
Acc = 0, Ay = 1 and /i = with N — 400 starting from (A^ and h 
kept fixed) (a) A^ = -0.2, ^ = 1, (b) A^ = 0.2, q^l. 



the odd revivals for even parity {q = 0). This can be under- 
stood easily by noting that for k = tt/N and t = ^ \ \^^^, 
we get Xk ^ ^ (as can be seen in the fidelity) and 2tAk ~ f , 
canceling the LE (eq. ([28])). 

Starting with a small non-zero value for A^^ and quenching 
the system to this critical point, the LE will behave roughly 
in the same way (Fig. [S]). However, we get two small peaks 
before the first revival. Strictly speaking, they will not be re- 
vivals according to the definition because they are less than 
two standard deviations away from the mean value of the LE. 
We can understand them as dynamical responses given by the 
Bogoliubov quasiparticles. 

We can also quench to the neighboring critical points along 
the critical curve A^ — A^^ A^y — 1 = 0, see (Fig. [9|). This means 
that both the starting and quenching Hamiltionian lie on this 
curve. In this case, the main difference will be the revival 
time, that will be exactly one third of the one found in the 
previous quenches. Note that we do not get this phenomenon 
for the critical line Xy = —Xx -\- 1, h = 0. 

In 1 16], the phenomenon of the revivals after a quantum 
quench is interpreted as a recombination of the fastest quasi- 
particles in the system. In general, the fastest excitations in 
the system have a speed that is upper bounded by the Lieb- 
Robinson speed ^lr |[40H42 | . This upper bound gives a lower 
bound to the revival time Tj-ev ^ Following (431, we find 

^LR — 3.2e/>/2 = 6.15 for the critical point A^ = 1, A^^ = 
and h = and therefore the maximum speed of quasi parti- 
cles given by the maximum group velocity 29/eA/e|max = 6 



(Fig. |9(d)| ) is compatible with the Lieb-Robinson bound. 
Conclusions. — In this paper, we studied the phase diagram 
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Figure 9. (Color online) LE for the quench to the critical point A^c = 
0, A^ = 1 and /i = along the critical line A^ — A^A^ — 1 = 
with N = 400 starting from (A^ along the critical line and /i = 
kept fixed ) (a) A^c = 0.1, ^ = 0. Notice that the revivals with odd 
parity are destroyed by interference (b) A^ = 0.1, ^ = 1. (c) Revival 
time for critical quenches to the same critical point starting in the 
neighboring ground states. Notice the critical line is detected by the 
revival time, (d) Group velocity 2dkAk for the critical point. The 
maximum value is 2^feAfc|max — 6. 



and quench behavior of the Cluster-XF model, a spin chain 
where the usual XY interactions in a transverse field are com- 
peting with a cluster three-body term. This model also de- 
scribes the effective behavior of the edge in a 2D fermionic 
symmetry protected topological state with Z2 symmetry |51 1. 
The Cluster-Xy model is exactly solvable by standard tech- 
niques, and the study has been conducted using the tools of 
fidelity susceptibility and Loschmidt Echo (LE). This model, 
inspired by proposed implementations of quantum compu- 
tation, provides a new benchmark with an interesting phe- 
nomenology and a much richer phase space coming from the 
competition of the different interactions. We were able to 
characterize the critical regions and the distribution of phases 
using the quantum geometric tensor. We found that the phase 
diagram is completely characterized by the fidelity. Notewor- 
thy, the ground state of some of the critical points present a 
large overlap with the ground state and few-excitations sub- 
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spaces of neighboring non-critical regions. The behavior 
away from equiUbrium is also non trivial. We showed that 
different critical points have qualitatively different effects on 
the LE. The long-time structure and the revival times of the 
LE depend on the initial phase of the quantum quench and the 
final critical point. This provides further phenomenology for 
the study of generic responses to critical quantum quenches. 

In [ 24, 44], it was shown that the Cluster-Ising model with 
open boundary conditions has a four fold degenerate ground 
space, which possesses symmetry protected topological or- 
der B31fl7ll . We expect that this model has similar features 
1341 , though the presence of the non trivial phases between 
the cluster phase and the ferromagnetic phase makes the situ- 
ation more complicated. A promising route to the characteri- 
zation of topological orders is the study of their entanglement 
spectrum [4^ ,49J . In a symmetry protected one-dimensional 
spin-one chain in the Haldane phase, the topological order is 



revealed in a double degeneracy of the entanglement spectrum 
1 50 1 . It would hence be interesting to study the entanglement 
spectrum of the cluster-XF model, to gain more insight into 
the properties of the symmetry protected topological order. In 
particular, it would be interesting to study the robustness of 
the information encoded in such ground space after a quan- 
tum quench, and whether the quench breaks or not the sym- 
metry that protects the topological order. Finally, it would 
be interesting to study the model in presence of disorder. In 
order to obtain reliable results, more sophisticated numerical 
techniques may be used such as matrix product states. 
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